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Abstract 

Medical imaging includes different modalities and processes to visualize the interior of human 
body for diagnostic and treatment purpose. However, one of the most common degradations in 
medical images is their poor contrast quality and noise. The existence of several objects and the 
close proximity of adjacent pixels values make the diagnostic process a daunting task. The idea of 
image enhancement techniques is to improve the quality of an image. In this study, morphological 
transform operation is carried out on medical images to enhance the contrast and quality. A disk 
shaped mask is used in Top-Hat and Bottom-Hat transform and this mask plays a vital role in the 
operation. Different types and sizes of medical images need different masks so that they can be 
successfully enhanced. The method shown in this study takes a mask of an arbitrary size and 
keeps changing its size until an optimum enhanced image is obtained from the transformation 
operation. The enhancement is achieved via an iterative exfoliation process. The results indicate 
that this method improves the contrast of medical images and can help with better diagnosis. 
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1. Introduction 

The field of digital image processing refers to processing digital images by means of a digital computer. In im- 
age processing operations, both the input and the output are images. Medical Image Processing and its applica- 
tions in Computer Assisted Diagnoses (CAD) and therapy (e.g. Computer Assisted Surgery — CAS) are of in- 
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creasing importance in modern medicine. Medical Image Enhancement is the art of examining images for iden- 
tifying objects and judging their significance. 

Most of the medical images (X-ray, CT, and MRI) have very low contrast and the challenge is to sharpen 
them. The effort on medical image enhancement has been focused mostly to improving visual perception of im- 
ages that are unclear because of blur. Edges are the representations of the discontinuities of image intensity 
functions. For processing these discontinuities in an image, a good edge enhancement technique is essential. 

Medical images enhancement is an important work for object recognition of the human organs such as lungs 
and ribs, and it is an essential pre-processing step in medical image segmentation [1] [2]. The work of contrast 
enhancement decides the result of the final processed image. Conventionally, edge is detected according to some 
early brought forward algorithms like Sobel algorithm, Prewitt algorithm and Laplacian of Gaussian operator [3], 
but in theory they belong to the high pass filtering, which are not fit for noise medical image edge detection be- 
cause noise and edge belong to the scope of high frequency. In real world applications, medical images contain 
object boundaries and object shadows and noise. Therefore, they may be difficult to distinguish the exact edge 
from noise or trivial geometric features. Mathematical morphology is a new mathematical theory which can be 
used to process and analyze the images [4] -[7]. It provides an alternative approach to image processing based on 
shape concept. 

Medical image enhancement technologies have attracted much attention since advanced medical equipments 
were put into use in the medical field. Enhanced medical images are desired by a surgeon to assist diagnosis and 
interpretation because medical image qualities are often deteriorated by noise and other data acquisition devices, 
illumination conditions, etc. Also targets of medical image enhancement are mainly to solve problems of low 
contrast and the high level noise of a medical image. Medical image enhancement technologies have attracted 
many studies [8]. The basic color space (RGB) conversion to other color spaces is important mainly concen- 
trated on the edges enhancement of medical color image. 

In the mathematical morphology theory, images are treated as sets, and morphological transformations which 
derived from Minkowski addition and subtraction are defined to extract features in images [9]. As the perfor- 
mance of classic edge detectors degrades with noise, morphological edge detector has been studied [10]. In this 
paper, a novel mathematical morphology edge detection algorithm is proposed based on using multi shape ap- 
proach whose information at various shapes is integrated to detect medical image edge. MATLAB was used 
here as the image processing tool like other research article [11]. 

2. Medical Image Enhancement 

Medical imaging is the technique and process of creating visual representations of the interior of a body for 
clinical analysis and medical intervention. Medical imaging seeks to reveal internal structures hidden by the skin 
and bones. It also makes possible the identification of abnormalities. 

Medical imaging encompasses different imaging modalities and processes to image the human body for di- 
agnostic and treatment purposes and therefore plays an important role in initiatives to improve public health for 
all population groups. It is crucial at all major levels of health care. In public health and preventive medicine as 
well as in both curative and palliative care, effective decisions depend on correct diagnoses. 

Image enhancement is a process principally focuses on processing an image in such a way that the processed 
image is more suitable than the original one for the specific application. The word “specific” has significance. It 
gives a clue that the results of such an operation are highly application dependent. In other words, an image en- 
hancement technique that works well for X-ray topographic images may not work well for MR images. 

Image enhancement approaches fall into two broad categories: spatial domain methods and frequency domain 
methods. The term spatial domain refers to the image plane itself, and approaches in this category are based on 
direct manipulation of pixels in an image. Frequency domain processing techniques are based on modifying the 
Fourier transform of an image. The process can be expressed as: g(x, y) = T[f (x, y)], where f(x, y) is the input 
image, g(x, y) is the processed image, and T is an operator on f defined over some neighborhood of (x, y). The 
principal approach in defining a neighborhood about a point (x, y) is to use a square or rectangular subimage 
area centered at (x, y). The center of the subimage is moved from pixel to pixel starting, say, at the top left cor- 
ner. The operator T is applied at each location (x, y) to yield the output, g, at that location. The process utilizes 
only the pixels in the area of the image spanned by the neighborhood. 

One method to enhance a medical image is wavelet transformation [12]. In this method, the high-frequency 
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sub-images are decomposed using the Haar wavelet transform. Then, noise in the high-frequency sub-bands is 
reduced using soft-thresholding. Finally, the enhanced image is obtained using inverse wavelet transformation. 

Another common method for medical image enhancement is histogram equalization [13] that enhances the 
contrast of image by increasing distribution of gray levels. This method does not necessarily obtain good results 
for all areas of an image because contrast enhancement may damage the image and the border areas. Because of 
this, there are different generalizations of this method to improve its performance [14] [15]. An algorithm for 
improving abdominal ultrasound images is proposed based on combination of histogram equalization and wave- 
let transformation in [16]. This algorithm improves edges and surroundings of abdominal walls and has real time 
performance in dynamic applications. Another method for medical image enhancement is Gamma correction 
[17]. In this method, Gamma values of individual pixels are locally optimized by minimizing the homogeneity 
of co-occurrence matrix of the original image. The Gamma correction method enhances dynamic range and im- 
proves the image. 

Morphological filters [18] are used for sharpening medical images. In this method, after locating edges by 
gradient-based operators, a class of morphological filter is applied to sharpen the existing edges. In fact, mor- 
phology operators, through increasing and decreasing colors in different parts of an image, have an important 
role in processing and detecting various existing objects in the image. Locating edges in an image using mor- 
phology gradient is an example that has comparable performance with that of classic edge-detectors such as 
Canny and Sobel [19]. 


3. Morphological Transform 

Methods of mathematical morphology act based on the structural properties of objects. These methods use ma- 
thematical principles and relationships between categories to extract the components of an image, which are 
useful in describing the shape of zones. Morphological operators are nonlinear, and two sets of data are their in- 
put. The first set contains the original image and the second one describes the structural element (mask). The 
original image is binary or in gray level and the mask is a matrix containing zero and one values. 

It is after applying the final image to the morphological operators that a new value for each pixel is obtained 
through sliding the mask on the original image. Value 1 in each mask indicates effectiveness and value 0 indi- 
cates ineffectiveness in the final image. Different formats can be selected to form a mask. Figure 1 shows a 
disk-shaped mask with radius of 4 (9 * 9 matrix). 


3.1. Morphological Operators 

If A (x, y) and B (u, v) describe the gray-level image matrix and the structural element matrix respectively, ero- 
sion and dilation operators are defined as (1) and (2): 

A0B = min |a(x + u,y + v)- B(u,v)| (1) 

A® B = max |A(x + u,y + v)-B(u,v)} (2) 

The erosion operator reduces the size of objects. This operator increases the size of holes in an image and re- 
moves very small details of that image. Removing bright areas under the mask makes the final image looks 
darker than the original image. The dilation operator acts in reverse; in other words, it increases and decreases 
the size of objects and holes in the image respectively. The opening operator is equivalent to the application of 
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Figure ] Disk-shaped structural element (mask) with radius of 4. 
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the erosion and dilation operations on the same image respectively (Equation (3)) while the closing operator acts 
in reverse (Equation (4)): 


A°B = (AGE)® B 

( 3 ) 

A-B = (A©B)©B 

( 4 ) 


The opening operator removes weak connections between objects and small details while the closing operator 
removes small holes and fills cracks. 

3.2. Selecting a Proper Mask 

Selecting a mask in proper shape and size to take morphological actions has a key role in achieving desired re- 
sults and reducing calculation time. In general, the shape and size of a mask are arbitrarily selected; however, 
the selected mask should be in appropriate shape and size for various diagnosis purposes. Disk-shaped masks 
(Figure 1) are more commonly used for medical images than other masks. As stated before, since disk-shaped 
masks are independent of changes in rotation, they are chosen for medical images. Since big or small masks 
strengthen or weaken various parts of an image, it is impossible to gather detailed information on the contrast of 
different images using only one structural element. This is why one mask in a particular shape and size may not 
appropriate for other applications [9]. In the proposed method, the change in shape and size of the mask contin- 
ues until an appropriate result obtained. It should be mentioned that past experiences have key roles in selecting 
proper masks to take morphological actions. 

3.3. Top-Hat Transforms 

These transforms are used to enhance the contrast of images through morphological methods and are in two 
general types: Top-Hat transform is obtained by subtracting the opening of the original image from the image 
itself (Equation (5)) and Bottom-Hat transform is obtained through subtracting the original image from its clos- 
ing (Equation (6)) [20]: 

Top-Hat (A) = A TH = A - ( A o B) (5) 

A q = A + ATH - ABH (6) 

Top-Hat and Bottom-Hat transforms are generally known as Open Top-Hat or White Top-Hat and Close Top- 
Hat or Black Top-Hat respectively. In many papers, Top-Hat is used to refer to both kinds of hat transforms. These 
transforms are named Top-Hat after Cylinder Hat shown in Figure 2. Figure 2(a) shows an image of the cross- 
section of the transform and Figure 2(b) describes defined parameters. In this figure, h indicates the difference in 
contrast and acts as a threshold, and d 0 and d z determine radius of the circles [21]. 

According to Equation (5), since the opening operator leaves a background of the image, it is expected that Top- 
Hat transform removes the image background. This transform acts like a high-pass filter and extracts the bright 
areas of the image (with contrast not less than h ) which are smaller than the mask. Bottom-Hat transform also re- 
moves the background of the image and leaves some dark areas of the image which are smaller than the mask 




(a) (b) 

Figure 2 . General idea in Top-Hat transform, (a): Cross-section of the transform; 
(b): Parameters of the transform. 
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itself [22]. It is possible to add the bright areas (the results of the opening operator) to the image and subtract the 
dark areas (the results of the closing operator) from it. As a result, there will be an enhancement in the contrast 
between bright and dark areas: 

A q = A + ATH - ABH (7) 


4. Methodology 

In the research Morphological Transform was implementing in MATLAB. At first we read the original image 
and calculate its CNR. We will require this CNR value to compare in the next steps. 

In the next step, we set an arbitrary structural element. We transform the original image and calculate CNR 
and PSNR values and store them. We compare this CNR value with the original image CNR value and if latter is 
larger we increase the size of structural element. Then we transform and calculate CNR and PSNR again and 
compare the current CNR value with the previous one. If the current value is larger we increase the size of 
structural element again. This goes on like a loop until we reach a point where the last CNR value is smaller 
than the previous one. We then set a structural element that is smaller than the last one and also set new end 
value and step of the loop to find out the structural element for whom the image has largest CNR value, meaning 
the image enhancement is optimum. The value we get at last is the size of our structural element for whom the 
original image is getting enhanced most. The number of steps and time taken to complete the operation is de- 
pendent on a few factors like size of the image, image source, initial condition, amount of noise etc. Figure 3 
shows the flowchart of the whole process. 

5. Results and Discussion 

From Figure 4(a) and Figure 4(b) we can see that optimum enhanced image is much clearer than the original 
image. This enhancement in contrast is the result of morphological transform carried out on the original image. 
We tried several masks and mask with a radius of 46 yielded the best result in this scenario. From Figure 4(c) 
and Figure 4(d) we can see that these have enhanced contrast than the original image but none of them are op- 
timum. Optimum enhanced image is clearer than both of them. Also, image enhanced with mask radius 70 has 
more darker spots along the edges although it worked better than mask with radius 30. 

From Figure 5 we can see that optimum enhanced image has more uniform intensity distribution than other 
images. It has its darker pixels optimized and distributed over the whole image as a result of morphological 
transform carried out on the original image. We also see that mask with a radius of 46 provided best results. 
Figure 6 represents the CNR values at different sizes of the structuring element. The highest CNR value is 
achieved when the structuring element has a radius of 46. 

Figure 7 demonstrated IVU bilateral elongated nephrograms with minimal irregularity of the renal outline. 
There is complete bilateral duplication of the renal collecting system to include duplication of the renal pelvis 
and uterus. The bilaterally duplicated uterus independently enters the bladder, and thus 4 UVJs are demonstrated. 
A small ureterocele of the right lower pole uterus is demonstrated at the confluence with the bladder. From 
above two images we can see that optimum enhanced image is much clearer than the original image. This en- 
hancement in contrast is the result of morphological transform carried out on the original image. We tried sever- 
al masks and mask with a radius of 16 worked best for this particular image. 

From Figure 7(c) and Figure 7(d) we can see that both images are clearer than the original and transformed 
image-2 shows better contrast than Transformed image-1 but has started growing darker pixels and hence it is 
not accepted as the final image. The morphological transform carried out on the original image increased the 
overall contrast and in the way also created some dark pixels which made portions of the second image opaque 
and hard to determine the actual structure. 

From Figure 8 we can see that optimum enhanced image has uniform intensity distribution than other images. 
It has its darker pixels optimized and distributed over the whole image as a result of morphological transform 
carried out on the original image. We also see that mask with a radius of 46 provided best results. There are less 
number of peaks indicating that overall contrast is increased (Table 1). 

We can see from Figure 9 that all three transformed images are clearer than the original but only one of them 
is optimum. Other two images have either less contrast than the optimum image or increased amount of dark 
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(a) (b) (c) (d) 


Figure 4. (a) Original image. This is an X-ray image of chest, abdomen and lower abdomen. Size = 2043 pixel 
by 2494 pixel CNR = 8.0996; (b) Optimum enhanced image where we used a structural element of radius 46. 
CNR = 16.2359, PSNR = 23.6968; (c) Transformed image-1. This image is not optimally enhanced. The 
structural element radius = 30, CNR = 15.5699, PSNR = 26.1075 and (d) Transformed image-2. This image is 
not optimally enhanced either when structural element radius = 70, CNR = 16.1918, PSNR = 22.6215. 



Figure 5. (a) Original image histogram; (b) Optimum enhanced image histogram. Structural element size = 46; 
(c) Transformed image-1 histogram. Structural element size = 30 and (d) Transformed image-2 histogram. 
Structural element size = 70. 




(a) (b) 

Figure 6. (a) CNR values at different sizes of the structuring element; (b) Highest CNR value is achieved when the structur- 
ing element has a radius of 46. 
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(a) (b) (c) (d) 


Figure 7. (a) Original image (X-ray of renal system). Size = 800 pixel by 1089 pixel CNR = 47.8057; (b) Opti- 
mum enhanced image. This image is optimally enhanced when we used a structural element of radius 16 CNR = 
51.2314, PSNR = 28.7180; (c) Transformed image-1. This image is not optimally enhanced. We used a structural 
element of radius 30. CNR = 48.3217, PSNR = 25.2594 and (d) Transformed image-2. This image is not opti- 
mally enhanced. Here structural element radius = 8, CNR = 48.5589, PSNR = 32.2399. 
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(c) (d) 

Figure 8. (a) Original image histogram; (c) Transformed image- 1 histogram. Structural element size = 30; (b) 
Optimum enhanced image histogram. Structural element size = 16; (d) Transformed image-2 histogram. Structur- 
al element size = 8. 


Table 1 Showing CNR and PSNR values after transforming the image with various sized structural elements. 


Structural element radius 

10 

12 

14 

16 

18 

20 

30 

CNR 

49.9110 

50.1925 

51.1017 

51.2314 

51.1416 

51.3295 

48.3271 

PSNR 

31.0268 

30.2598 

29.3051 

28.7180 

28.1790 

27.4393 

25.2595 
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spots which make it harder to see the internal structure. This enhancement in contrast is the result of morphological 
transform carried out on the original image. We tried several masks and mask with a radius of 36 worked best. 

From Figure 10 we can see that optimum enhanced image has uniform intensity distribution than other im- 
ages. It has fewer amounts of dark pixels and increased number of bright pixels. This enhancement in contrast is 
the result of morphological transform carried out on the original image. We tried several masks and mask with a 
radius of 36 yielded the best result (Table 2). 

From Figure 11 we can see that all three processed images better contrast than the original image but opti- 
mum enhanced image has its edges preserved which is important. Transformed image-2 is also enhanced but it 
has dark pixels along the edges of the left ear which has blurred its edge (Table 3 and Table 4). 

From Figure 12 we can see that optimum enhanced image has uniform intensity distribution than other im- 
ages. It has its darker pixels optimized and distributed over the whole image as a result of morphological trans- 
form carried out on the original image. There is less number of peaks indicating that overall contrast is increased. 
We tried several masks and mask with a radius of 154 worked best in this scenario. 



Figure 9. (a) Original image (MRI image of brain) showing Langerhans cell histiocytosis. Size = 512 pixel by 512 pixel 
CNR = 57.391; (b) Optimum enhanced image. This image is optimally enhanced. Structural element radius = 36, CNR = 
79.232, PSNR = 25.5522; (c) Transformed image-1. This image is not optimally enhanced. Structural element radius = 10, 
CNR = 76.9586, PSNR = 28.6038 and (d) Transformed image-2. This image is not optimally enhanced. Structural element 
radius = 70, CNR = 77.555, PSNR = 23.6408. 





(c) 


(d) 


Figure 10. (a) Original image histogram; (b) Optimum enhanced image histogram. Structural element size = 36; (c) Trans- 
formed image- 1 histogram. Structural element size = 10 and (d) Transformed image-2 histogram. Structural element size = 
70. 
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(a) (b) (c) (d) 


Figure 11. (a) Original image (X-ray) showing Langerhans cell histiocytosis. The bones are well-mineralized. No radio- 
graphic evidence of rickets is seen. No classic metaphyseal lesion is seen. Size = 1566 pixel by 1602 pixel CNR = 20.4178; 
(b) Optimum enhanced image. This image is optimally enhanced. Structural element radius = 154, CNR = 41.2622, PSNR = 
20.0363; (c) Transformed image-1. This image is not optimally enhanced. Structural element radius = 120, CNR = 41.1028, 
PSNR = 21.0572 and (d)-Transformed image-2. This image is not optimally enhanced. Structural element radius = 190, CNR 
= 41.056, PSNR = 18.5121. 
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Figure 12 (a) Original image histogram; (b) Optimum enhanced image histogram. Structural element size = 154; (c) Trans- 
formed image-1 histogram. Structural element size = 120 and (d) Transformed image-2 histogram. Structural element size = 
190. 


Table 2 Showing CNR and PSNR values after transforming the image with various sized structural elements. 


Structural element radius 

10 

20 

30 

32 

34 

36 

38 

40 

50 

CNR 

76.9586 

78.6143 

79.0811 

79.1315 

79.1954 

79.2320 

79.2050 

79.1770 

79.0347 

PSNR 

28.6039 

26.7137 

25.8554 

25.7543 

25.6152 

25.5522 

25.4643 

25.3485 

24.7570 
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Table 3 CNR and PSNR values after transforming the image with various sized structural elements (10 to 100). 


Structural Element Radius 

10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

CNR 

36.8571 

38.3132 

39.1797 

39.6674 

39.9473 

40.1759 

40.4481 

40.6399 

40.802 

40.944 

PSNR 

31.9395 

27.9438 

25.3337 

24.2909 

23.7417 

23.4342 

22.9278 

22.3743 

22.024 

21.695 


Table 4 CNR and PSNR values after transforming the image with various sized structural elements (110 to 170). 


Structural Element Radius 

110 

120 

130 

140 

150 

152 

154 

156 

160 

170 

CNR 

41.0495 

41.1028 

41.1645 

41.1915 

41.2325 

41.2561 

41.2622 

41.2437 

41.232 

41.198 

PSNR 

21.2969 

21.0572 

20.6421 

20.4126 

20.1689 

20.0898 

20.0363 

19.8874 

19.768 

19.101 


6. Conclusion 

We observed the application of morphological transform in image manipulation. Contrast of various medical 
images can be enhanced using this method. We implemented this technique using Image Processing Toolbox in 
MATLAB. The processed images are created by adding the Top-Hat transformation of the original image to the 
original image and then subtracting the Bottom-Hat transformation from it. For both these transforms, we used a 
structural element of arbitrary size at first and then increased it stepwise to achieve a value that gave optimum 
contrast enhancement for the image that we were working with. We also calculated contrast-to-noise ratio to 
compare two processed images. The newly created images show better contrast as we see than the original im- 
ages. This is due to the morphological operations carried out on the pixels. This technique can be used to en- 
hance various medical images that can help medical practitioners with better diagnosis and treatment which are 
invaluable to improve public health. 
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